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1 Introduction 

Evidence for dark matter in the universe is available from 
a wide range of observational data. In addition to the clas- 
sic evidence from galactic rotation curves [T] , the analysis 
of the cosmic microwave background anisotropics leads to 
the conclusion that the curvature of the universe is close 
to zero indicating that the sum of the fractions of critical 
density, fi, in matter and a cosmological constant (or dark 
energy) is close to one [2] . When combined with a variety 
of data including results from the analysis of type la su- 
pernovae observations [3] and baryon acoustic oscillations 
[3] one is led to the concordance model where f2rn ~ 0.23 
and f^A ~ 0.73 with the remainder (leading to f2tot = 1) 
in baryonic matter. Other dramatic pieces of evidence can 
be found in combinations of X-ray observations and weak 
lensing showing the superposition of dark matter (from 
lensing) and ordinary matter from X-ray gas and from 
the separation of baryonic and dark matter after the col- 
lision of two galaxies as seen in the Bullet cluster [5]. For 
a more complete discussion see [7]. 

Here, I will adopt the results of the three-year data 
from WMAP |2j which has determined many cosmologi- 
cal parameters to unprecedented precision. Of particular 
interest is the determination of the total matter density 
(relative to the critical density), fimh^ and the baryonic 
density, Slbh^. In the context of the ylCDM model, the 
WMAP only results indicate 



i7„/i^ = 0.1265; 



-0.0081 
-0.0080 



= 0.0223 ± 0.0007 (1) 



The difference corresponds to the requisite dark matter 
density 



^CDMh^ = 0.1042~''n'nnfin 



0.0080 



(2) 



or a 2cr range of 0.0882 - 0.1204 for QcDAih'^ 



^ To be published in "Supersymmetry on the Eve of the 
LHC" a special volume of European Physical Journal C, Par- 
ticles and Fields (EPJC) in memory of Julius Wess. 



2 Constrained IVISSM models 

To construct the supersymmetric standard model ^ we 
start with the complete set of chiral fermions needed in the 
standard model, and add a scalar superpartner to each 
Weyl fcrmion so that each field in the standard model 
corresponds to a chiral multiplet. Similarly we must add 
a gaugino for each of the gauge bosons in the standard 
model making up the gauge multiplets. The minimal su- 
persymmetric standard model (MSSM) [S] is defined by 
its minimal field content (which accounts for the known 
standard model fields) and minimal superpotential nec- 
essary to account for the known Yukawa mass terms. As 
such we define the MSSM by the superpotential 

W = e„ [veHiVe^ + ydH{Q'd- + VuHlQ^u^] + (3) 



where 



(4) 



In ([3]), the indices, {ij}, are SU(2)i doublet indices. The 
Yukawa couplings, y, are all 3 x 3 matrices in genera- 
tion space. Note that there is no generation index for the 
Higgs multiplets. Color and generation indices have been 
suppressed in the above expression. There are two Higgs 
doublets in the MSSM. This is a necessary addition to 
the standard model which can be seen as arising from the 
holomorphic property of the superpotential. That is, there 
would be no way to account for all of the Yukawa terms 
for both up-type and down-type multiplets with a single 
Higgs doublet. To avoid a massless Higgs state, a mixing 
term must be added to the superpotential. 

In order to preserve the hierarchy between the elec- 
troweak and GUT or Planck scales, it is necessary that the 
explicit breaking of supersymmetry be done softly, i.e., by 
the insertion of weak scale mass terms in the Lagrangian. 
This ensures that the theory remain free of quadratic di- 
vergences [To]. The possible forms for such terms are 



^so 



■ft 
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^{Ayf''^,(l)j^k + h.c. (5) 



where the are gaugino masses, w? are soft scalar 
masses, B is a bilinear mass term, and A is a trilinear 
mass term. Masses for the gauge bosons are of course for- 
bidden by gauge invariance and masses for chiral fermions 
are redundant as such terms are explicitly present in M*^ 
already. For a more complete discussion of supersymemtry 
and the construction of the MSSM see ITII. 



2.1 Neutralinos and the Relic Density 

There are four neutralinos, each of which is a linear combi- 
nation of the i? = — 1 neutral fermions [12]: the wino , 
the partner of the 3rd component of the SU{2)l gauge 
boson; the bino, B] and the two neutral Higgsinos, Hi 
and H2 ■ The mass and composition of the LSP are deter- 
mined by the gaugino masses, ^, and tan/3. In general, 
neutralinos can be expressed as a linear combination 



(6) 



The solution for the coefficients a,/3,7 and 6 for neutrali- 
nos that make up the LSP can be found by diagonalizing 
the mass matrix 



/ 



{W\B,H"i,H^) 



M2 


-32fl 

1 g2V2 

\ V2 





A/i 

-giv2 
V2 



V2 





g2V2 

-giV2 
V2 





B 
Hi 

\H"2J 



(7) 

where Mi{M2) is a soft supersymmetry breaking term giv- 
ing mass to the U(l) (SU(2)) gaugino(s). 

The relic abundance of LSP's is determined by solving 
the Boltzmann equation for the LSP number density in 
an expanding Universe, 



— = {(Tv)(n -no) 



(8) 



where rig is the equilibrium number density of neutralinos. 
By defining the quantity / = n/T^, we can rewrite this 
equation in terms oi x = T jm^, as 



dx 



-^Gnn) if - f) (9) 



of interest, the LSP is a bino and the annihilation pro- 
ceeds mainly through sfermion exchange. Because of the 
p-wave suppression associated with Majorana fermions, 
the s-wave part of the annihilation cross-section is sup- 
pressed by the outgoing fermion masses. This means that 
it is necessary to expand the cross-section to include p- 
wave corrections which can be expressed as a term pro- 
portional to the temperature if neutralinos are in equilib- 
rium. Unless the neutralino mass happens to lie near near 
a pole, such as ~ mz/2 or mh/2, in which case there 
are large contributions to the annihilation through direct 
s-channel resonance exchange, the dominant contribution 
to the BB annihilation cross section comes from crossed 
i-channel sfermion exchange. 

Annihilations in the early Universe continue until the 
annihilation rate F ~ avn~f^ drops below the expansion 
rate. The final neutralino relic density expressed as a frac- 
tion of the critical energy density can be written as [12] 



Q^h'' ~ 1.9 X 10" 



GeV 



axf + \bx'\ 



(10) 



where [T^/T^Y accounts for the subsequent reheating of 
the photon temperature with respect to x, due to the an- 
nihilations of particles with mass m < XfUi^ |15j and 
Xf = Tf / is proportional to the freeze-out tempera- 
ture. The coefficients a and b are related to the partial 
wave expansion of the cross-section, av = a + bx + . . .. 
Eq. (|10l ) results in a very good approximation to the relic 
density expect near s-channel annihilation poles, thresh- 
olds and in regions where the LSP is nearly degenerate 
with the next lightest supersymmetric particle [16]. 

When there are several particle species «, which are 
nearly degenerate in mass, co-annihilations are important. 
In this case[l6j, the rate equation ([8]) still applies, provided 
n is interpreted as the total number density. 



(11) 



riQ as the total equilibrium number density, 

no = ^no,i , (12) 



and the effective annihilation cross section as 



(CcffWrcl) 



E 



no,inoj 



(13) 



The solution to this equation at late times (small x) yields 
a constant value of /, so that n (x T^. 

The technique (13] used to determine the relic density 
is similar to that for computing the relic abundance of 
massive neutrinos [14] with the appropriate substitution of 
the cross section. The relic density depends on additional 
parameters in the MSSM beyond Mi,M2,/u, and tan/3. 
These include the sfermion masses, and the Higgs 
pseudo-scalar mass, niA- To determine the relic density 
it is necessary to obtain the general annihilation cross- 
section for neutralinos. In much of the parameter space 



In eq. ([9]), is now understood to be the mass of the 
lightest sparticle under consideration. 



2.2 The CMSSM 

In its generality, the minimal supersymmetric standard 
model (MSSM) has over 100 undetermined parameters. 
But in addition to relieving the helplessness of an anal- 
ysis with so many free parameters, there are good argu- 
ments based on grand unification [17j and supergravity 
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|18| which lead to a strong reduction in the number of 
parameters. I will assume several unification conditions 
placed on the supersymmetric parameters. In all models 
considered, the gaugino masses are assumed to be unified 
at the GUT scale with value, TOi/2, as are the trilinear 
couplings with value Aq. Also common to all models con- 
sidered here is the unification of all soft scalar masses set 
equal to toq at the GUT scale. With this set of boundary 
conditions at the GUT scale, we can use the the radiative 
electroweak symmetry breaking conditions by specifying 
the ratio of the two Higgs vacuum expectation values, 
tan P, and the mass, Mz, to predict the values of the Higgs 
mixing mass parameter, /x and the bilinear coupling, B. 
The sign of fj, remains free. This class of models is often 
referred to as the constrained MSSM (CMSSM) pl,[20 | l2 H 
[21123]. In the CMSSM, the solutions for ^ generally lead 
to a lightest neutralino which is very nearly a pure B. 

Having set the boundary conditions at the GUT scale, 
renormalization group equations are run down to the elec- 
troweakscale. For example, at 1-loop, the RGEs for the 
gaugino masses are: 



dt 



(14) 



Assuming a common gaugino mass, mi/2 at the GUT scale 
as was discussed earlier, these equations are easily solved 
in terms of the fine structure constants. 



ai{t) 



ai{MGUT) 



This implies that 



Ml Mo 



9 9 

91 .92 



9l 



(15) 



(16) 



(Actually, in a GUT, one must modify the relation due to 
the difference between the U(l) factors in the GUT and 
the standard model, so that we have Mi = |^Af2.) At 
two loops, these relations are modified. 

In the figure below, an example of the running of the 
mass parameters in the CMSSM is shown. Here, we have 
chosen mi/2 = 250 GeV, mo = 100 GeV, tan/3 = 3, 
^0 — 0; and /i < 0. Indeed, it is rather amazing that 
from so few input parameters, all of the masses of the su- 
persymmetric particles can be determined. The character- 
istic features that one sees in the figure, are for example, 
that the colored sparticles are typically the heaviest in the 
spectrum. This is due to the large positive correction to 
the masses due to in the RGB's. Also, one finds that 
the B, is typically the lightest sparticle. But most impor- 
tantly, notice that one of the Higgs mass^, goes negative 
triggering electroweak symmetry breaking [24^ . (The neg- 
ative sign in the figure refers to the sign of the mass^, even 
though it is the mass of the sparticles which is depicted.) 

Although the CMSSM is often confused in name with 
mSUGRA, i.e. models based on minimal supergravity, 
[18], the latter employ two additional constraints [26]. In 
the simplest version of the t heor v |?7[[^[TS] where super- 
symmetry is broken in a hidden sector, the universal tri- 
linear soft supersymmetry-breaking terms are A^ — {i — 
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Fig. 1. RG evolution of the mass parameters in the CMSSM. 
I thank Toby Falk for providing this figure. 



\/3)too and bilinear soft supersymmetry-breaking term is 
So = (2 — V3)mo, i.e., a special case of a general relation 
between B and A, Bq = Aq — wlq. 

Given a relation between Bq and we can no longer 
use the standard CMSSM boundary conditions, in which 
"T'i/2> 'TfiQ, Aq, tan/3, and sgn{^) are input at the GUT 
scale with fi and B determined by the electroweak sym- 
metry breaking condition. Now, one is forced to input Bq 
and instead tan P is calculated from the minimization of 
the Higgs potential [26]. 

In addition, there is a relation between the gravitino 
mass and soft scalar masses, 111^/2 — ttlq. When elec- 
troweak symmetry breaking boundary conditions are ap- 
plied, this theory contains only mi/ 2, tiq, and Aq in addi- 
tion to the sign of ^, as free parameters. The magnitude 
of as well as tan (3 are predicted. 

In contrast to the mSUGRA model which contains one 
fewer parameter than the CMSSM, the remaining mod- 
els to be discussed below are less constrained than the 
CMSSM, in that they involve more free parameters. One 
is a GUT- less model [25" . It is GUT-less in the sense that 
the supersymmetry breaking masses are unified not at the 
GUT scale, but at an input scale chosen to be below the 
GUT scale. This input parameter, Mi„, is the one ad- 
ditional parameter relative to the CMSSM in the GUT- 
less models. Gauge coupling unification still occurs at the 
GUT scale. In the second class of models, the NUHM, the 
Higgs soft masses are not unified at the GUT scale [521 
[30] , In this class of models, both and the Higgs pseudo 
scalar mass become free parameters. 

In all of the models discussed below, I will assume un- 
broken R — parity and and hence the lightest supersym- 
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metric particle is stable. This will often, but not always 
be the neutralino 12 1 . 



3 The CMSSM after WMAP 

For a given value of tan/3, Aq, and sgn{ii), the result- 
ing regions of acceptable relic density and which satisfy 
the phenomenological constraints can be displayed on the 
mi/2 — mo plane. In Fig. [2k, the light shaded region corre- 
sponds to that portion of the CMSSM plane with tan (3 = 
10, Aq — 0, and /i > such that the computed relic den- 
sity yields the WMAP value given in eq. Q gS]. The bulk 
region at relatively low values of TOi/2 and mo, tapers off 
as my 2 is increased. At higher values of mo, annihilation 
cross sections are too small to maintain an acceptable relic 
density and f2^h^ is too large. Although sfermion masses 
are also enhanced at large mi/2 (due to RGE running), 
co-annihilation processes between the LSP and the next 
lightest sparticle (in this case the r) enhance the anni- 
hilation cross section and reduce the relic density. This 
occurs when the LSP and NLSP are nearly degenerate 
in mass. The dark shaded region has mf < and is ex- 
cluded. The effect of coannihilations is to create an allowed 
band about 25-50 GeV wide in mo for mi/2 ^ 950 GeV, 
or mi/2 ^ 400 GeV, which tracks above the mf ^ = 
contour [31 . 

Also shown in Fig. [2^ are the relevant phenomeno- 
logical constraints. These include the LEP limits on the 



tan p = 10 , p. > 



chargino mass: m^± > 104 GeV 32], on the selectron 
mass: mg > 99 GeV [33] and on the Higgs mass: m^ > 
114 GeV f31|. The former two constrain mi/2 and mo di- 
rectly via the sparticle masses, and the latter indirectly via 
the sensitivity of radiative corrections to the Higgs mass 
to the sparticle masses, principally m^^. FeynHiggs [35j 
is used for the calculation of m/j. The Higgs limit imposes 
important constraints principally on mi/2 particularly at 
low tan /3. Another constraint is the requirement that the 
branching ratio for 6 ^ 57 is consistent with the ex- 
perimental measurements f36^ . These measurements agree 
with the Standard Model, and therefore provide bounds 
on MSSM particles [37], such as the chargino and charged 
Higgs masses, in particular. Typically, the 6 — *■ 57 con- 
straint is more important for /i < 0, but it is also relevant 
for /i > 0, particularly when tan/3 is large. The constraint 
imposed by measurements of 6 — *■ 57 also excludes small 
values of mi/2- Finally, there are regions of the (mi/2, mo) 
plane that are favoured by the BNL measurement ^38] of 
(7^ — 2 at the 2-a level, corresponding to a deviation from 
the Standard Model calculation [39] using e+e~ data. 

Another mechanism for extending the allowed regions 
in the CMSSM to large is rapid annihilation via a 
direct- channel pole when m^ ~ ^m^i [T9l[2T] . Since the 
heavy scalar and pseudoscalar Higgs masses decrease as 
tan P increases, eventually 2m^ ~ ttia yielding a 'funnel' 
extending to large rrii /2 and mo at large tan /?, as seen 
in Fig. [2b- As one can see, the impact of the Higgs mass 
constraint is reduced (relative to the case with tan/3 — 10) 
while that of 6 — > 37 is enhanced. 
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Fig. 2. The (mi/2, nio) planes for (a) tan /3 = 10 and fi > 0, 

assuming Aq = 0, nit ~ 175 GeV and m,b{mb)'sM = 4.25 GeV. 
The near-vertical (red) dot-dashed lines are the contours ruh ~ 
114 GeV, and the near-vertical (black) dashed line is the con- 
tour m^± — 104 GeV. Also shown by the dot-dashed curve 
in the lower left is the corner excluded by the LEP bound of 
mg > 99 GeV. The medium (dark green) shaded region is ex- 
cluded by b sj, and the light (turquoise) shaded area is the 
cosmologically preferred region. In the dark (brick red) shaded 
region, the LSP is the charged fi. The region allowed by the 
E821 measurement of at the 2-a level, is shaded (pink) and 
bounded by solid black lines, with dashed lines indicating the 
l-o ranges. In (b), tan /? = 50. 



Shown in Fig.[3]are the WMAP lines [H] of the (mi/2, mo) 
plane for /i > and values of tan /3 from 5 to 55, in steps 
Z\(tan/3) = 5. We notice immediately that the strips are 
considerably narrower than the spacing between them, 
though any intermediate point in the (mi/2,mo) plane 
would be compatible with some intermediate value of tan /3. 
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The right (left) ends of the strips correspond to the maxi- 
mal (minimal) allowed values of mi/2 and hence m^. The 
lower bounds on TOi/2 are due to the Higgs mass con- 
straint for tan/3 < 23, but are determined by the 6 — > 57 
constraint for higher values of tan /?. 
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Fig. 3. The strips display the regions of the (mi/2, mo) 
plane that are compatible with the WMAP determination of 
Q^h^ and the laboratory constraints for fi > and tan/3 = 
5,10,15,20,25,30,35,40,45,50,55. The parts of the strips 
compatible with g^^ — 2 at the 2-cr level have darker shaa 



Finally, there is one additional region of acceptable 
relic density known as the focus-point region [40], which 
is found at very high values of toq. An example showing 
this region is found in Fig.[4l plotted for tan/3 = 10, > 0, 
and m* = 175 TeV. As mo is increased, the solution for 
at low energies as determined by the electroweak symme- 
try breaking conditions eventually begins to drop. When 
jjL < mi/2, the composition of the LSP gains a strong Hig- 
gsino component and as such the relic density begins to 
drop precipitously. As mo is increased further, there are 
no longer any solutions for /i. This occurs in the shaded 
region in the upper left corner of Fig. 3] The position of 
the focus point strip is very sensitive to the value of rrit 

m. 



4 Direct detection 

Direct detection techniques rely on the neutralino nucleon 
scattering cross-section. In general, there are two contribu- 
tions to the low-energy effective four-fermion Lagrangian 
which are not velocity dependent. These can be expressed 
as spin-dependent and scalar or spin-independent interac- 
tions. 
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mm (GeV) 

Fig. 4. As in Fig. Oi, where the range in mo is extended to 5 
TeV. In the shaded region at very high mo, there are no solu- 
tions for fj, which respect the low energy electroweak symmetry 
breaking conditions. 



which is to be summed over the quark flavours q, and the 
subscript i labels up- type quarks (i — 1) and down- type 
quarks (i = 2). The detailed expressions for a2i and a^i 
can be found in ■ 

The scalar cross section for a neutralino scattering on 
a nucleus with atomic number A and charge Z, can be 
written in terms of asi, 



where 



4m o 
^3 = — [ZU + {A- Z)U , 

TT 

ir is the reduced LSP mass, 
fp_ _ \^ Ap) I 2 (p) asg 

m _ 



q—u,d.s 



(18) 



(19) 



c,b,t 



(20) 



the parameters fji^ are defined by 

mpf^g = (p\mqqq\p) =mqBg, 

Itg = 1 ~ Eg=«,d,s frl 03], and /„ has a similar expres- 
sion. The needed matrix elements are determined in part 
by the 7r-nucleon S term, which is given by 



(21) 



and the strangeness contribution to the proton mass, 
2Bs 



y 



Bu + Bd 



1 



CTo 



(22) 



where ctq is the change in the nucleon mass due to the non- 
zero u, d quark masses, which is estimated on the basis of 
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Table 1. Matrix elements in terms of proton strangeness 



E (MeV) 


y 








45 


0.2 


0.020 


0.026 


0.117 


64 


0.44 


0.027 


0.039 


0.363 


36 





0.016 


0.020 


0.0 



octet baryon mass differences to be co = 36 ± 7 MeV [44] . 
Values of the matrix elements are given in the table for 3 
choices of E and hence y. As one can see, we expect there 
will be a strong dependence of the cross section on the 
proton strangeness, y. 

The spin-dependent (SD) part of the elastic x-nucleus 
cross section can be written as Q 

32 

TSD = —GlmlA^J{J + 1) , (23) 

TT 

where is again the reduced neutralino mass, J is the 
spin of the nucleus, 

^= l(ap(5p)+a„(5„)) , (24) 

and 

The factors parametrize the quark spin content of 

the nucleon and are only significant for the light (u,d,s) 
quarks. A combination of experimental and theoretical re- 
sults tightly constrain the linear combinations 

EE Z\(P) - Z\|iP^ = 1.2695 ± 0.0029 (26) 

and 

= Z\(,P) + Z\^P^ - 2zi(P) = 0.585 ± 0.025. (27) 

However, the individual Aq^^ are relatively poorly con- 
strained; using the recent COMPASS result 48J, 

Ai^ls = -0.09±0.01 (stat.)±0.02 (syst.) w -0.09±0.03 , 

(28) 

where we have conservatively combined the statistical and 
systematic uncertainties, we may express A^f^J as follows 
in terms of known quantities [49j : 

A[P^ = i (4"^ + "s"^) + ^iF^ = 0-84 ± 0.03 (29) 

A^f = i (4?) _ 4?)) + ^p) ^ -0.43 ± 0.03 (30) 

^ As with the SI cross section, this expression applies in the 
zero momentum transfer limit and requires an additional form 
factor for finite momentum transfer. This form factor may have 
a small but non-zero dependence on Op and an. 



The above two uncertainties and that of As , however, 
are correlated and one should instead use the indepen- 
dent quantities a^^"^ and Og'''' when appropriate. The pro- 
ton and neutron scalar matrix elements are related by an 
interchange of A^ and A^, or 

A(n)^^(p)_ /iW^^(p)^ and Z\(")=Z\(P). (31) 

In addition to the uncertainties stemming from St^n , it 
is important to bear in mind the astrophysical uncertain- 
ties inherent in direct detection experiments. Direct and 
indirect detection signals are proportional to the product 
of the local dark matter density po and the elastic cross 
section. Since the local dark matter density has not been 
directly measured, it is typically inferred from galactic dy- 
namics and iV-body simulations. By convention, experi- 
mental results are often presented for a neutralino density 
of 0.3 GeV/cm^, with the implicit understanding that the 
following significant uncertainties exist in this value. In 
the case of a smooth distribution of galactic dark mat- 
ter, po is estimated to be 0.2-0.4 GeV/cm^ [SOllST] for 
a spherical halo, but it may be somewhat higher, up to 
0.7 GeV/cm3 [S113], for an eUiptical halo; see Ref. [FT] 
for a discussion of the difficulties in determining this value. 
Models of the galaxy based upon hierarchical formation [54] , 
which do not assume a strictly smooth distribution of dark 
matter as do the above estimates, suggest that the lo- 
cal density may be as low as 0.04 GeV/cm'^, but that 
0.2 GeV/cm'^ is a more reasonable lower limit. One should 
also allow for the possibility that there may be some addi- 
tional source of cold dark matter, such as axions, in which 
case the neutralino density would be less than po- 

As noted earlier, regions in the ■m-i^/2,mQ plane with 
fixed tan/3 for which the relic density takes on values 
compatible with the WMAP determined cold dark matter 
density, form thin strips due to co-annihilations or rapid 
s-channel annihilations. As a consequence, for fixed tan/?, 
we can adjust mp as a function of mi/2 and calculate the 
elastic scattering cross section as a function of fTii/2. Fig. 
[5] [4£| shows the cross sections along the WMAP-allowed 
strips of the (mi/2,TOo) planes for tan/3 = 10 and 50 [22] 
[^[55] . Shown are both the spin dependent cross section 
for x~P and x ~ scattering (upper curves) and the spin 
independent cross section (lower curves, which are nearly 
superimposed at this scale). 

The SI cross sections for x-proton (cr;^p,si) and x-neutron 
((T^ii^i) scattering are typically very close, within ~5% of 
each other; close enough that they are virtually indistin- 
guishable in Fig.[5j Likewise, the spin-dependent x-proton 
("■xp.Sd) and x-neutron (cr^n.SD) scattering cross sections 
are similar, differing by at most a factor of 2 or 3. In gen- 
eral, the SD x-nucleon cross section is much larger than 
the SI one, by 0(10^ — 10"^) or more. However, we recall 
that the SI cross section for scattering off a nucleus, Eq. 
1181 contains a factor of the number of nucleons squared, 
whereas the SD cross section, Eq.[231 is proportional to the 
square of the spin, which does not grow with increasing 
nuclear mass. Consequently, heavy elements such as Ge 
and Xe are actually more sensitive to SI couplings than to 
SD couplings. 
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Fig. 5. The neutralino-nucleon scattering cross sections along the WMAP-allowed coanmhilation strip for tan/3 = 10 (left 
panel) and coannihilation/ funnel strip for tan/3 = 50 (right panel) 




Fig. 6. The spin-independent neutralino-nucleon scattering cross section ratio along the WMAP-allowed strips for tan /? = 10 
(left panel) and tan/3 — 50 (right panel) for several values of Sttn- Experimental constraints from CDMS J50j and XENONIO 
\57^ are also shown. 



In Fig. [HI the SI cross sections are shown along the 
WMAP allowed coannihilation strip for tan/3 = 10 and 
coannihilation/funnel strip for tan/3 — 50 for the St^m 
reference values of 36 MeV (no strange scalar contribu- 
tion), 64 MeV (central value), and 80 MeV {2-a upper 
bound) 39]. The CDMS and XENONIO hmits 
are also given. A factor of ~10 variation occurs in ct^at.si 
over these X'^jv reference values for any given model along 
the WMAP strip. 

Such large variations present difficulties in using any 
upper limit or possible future precision measurement of 
'^xN,si from a direct detection signal to constrain the 
CMSSM parameters. The present CDMS and XENONIO 
upper limits have (almost) no impact on the WMAP strip 
for tan/3 = 10(50), if one makes the very conservative as- 
sumption that (To — 36 MeV (y = 0). On the other hand, 
mi/2 ^ 200 GeV would be excluded for tan/3 = 10 if 



^■KN = 64 or 80 MeV. This excluded region would ex- 
tend to mi/2 ^ 300 GeV for tan/3 = 50 if S-^m = 64 
or 80 MeV. Thus, the experimental uncertainty in St^^ 
is already impinging on the ability of the present CDMS 
and XENONIO resuhs to constrain the CMSSM parame- 
ter space [49] . 



Looking to the future, a conjectural future measure- 
ment of CT^p.si = 4 X 10~^ pb would only constrain toi/2 
to the range 600 GeV < m^/^ < 925 GeV if tan^ = 10 
and 1100 GeV < m^/^ < 1400 GeV if tan/3 = 50, for the 
l-tJ i^TOv range of 64 ± 8 MeV. If smaller values of Z"™ 
are also considered, down to o-q = 36 MeV (y — 0), these 
constraints would weaken to 350 GeV < mi/2 < 925 GeV 
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and 550 GeV < mi/2 < 1400 GeV for tan/3 = 10 and 50, 
respectively, ruling out only the smallest values of mi/2 0- 
In Figure El we show the S-„m dependence of cr^Af,si 
[491158] for three benchmark models, which are essentially 
the well-studied benchmark models C, L, and M [59 . All 
of the benchmark models lead to relic densities within 
the WMAP preferred range of Qh"^ = 0.088 - 0.120 ;2 , 
and satisfy most phenomenological constraints [fj. Points 
C and L are points along the coannihilation strip (at low 
tan/3 = 10 and high tan/? = 50, respectively) where the 
masses of the neutralino and stau are nearly degenerate. 
Point M is in the funnel region where the mass of the 
neutralino is roughly half that of the Higgs pseudoscalar. 
From the minimal value for (co — 36 MeV) to the 
2-cr upper bound (80 MeV), <JxN,si varies by more than 
a factor of 10 (as much as a factor of 35 for model L). 
At these benchmarks and in other models of interest, for 
larger values of Z'^n, the majority of the contribution to 
/p in Eq. [19] comes from the strange quark term, with 

f^^ (xy (X S^t, - (To, so that ctxp.Si ^ i^™ - ctq)"^. Thus, 
the SI cross sections are particularly sensitive not just to 
Z'ttn and (To J but to their difference. For smaller values 
~ (To), the strange contribution no longer 
dominates, but a strong dependence of cf^n.si on does 
remain. 
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(MeV) 

Fig. 7. The spin-independent neutralino-nucleon scattering 
cross section as a function of X'^n for benchmark models C, 
L, and M. Note that (T^p.si and (Txn.si are nearly indistinguish- 
able at the scale used in this plot. 



In Fig. [8] we display the expected ranges of the spin- 
independent cross sections in the CMSSM when we sample 
randomly tan /3 as well as the other CMSSM parameters 



^ Moreover, as noted previously, however, detection signals 
only measure poa^N, so can only be determined from a 
signal to the precision that the local dark matter density is 
known. 

^ The exception being a possible failure to account for the 
discrepancy between the theoretical and experimental results 
for anomalous magnetic moment of the muon [381139) . Points 
C and L are consistent with {g — 2)^, whilst the contribution 
from point M is too small. 



[58] . All points shown satisfy all phenomenological con- 
straints and have f^x^^ ^^^^ than the WMAP upper limit. 
Since models with low fi^^^ '-■^'^ '^'^^ be excluded if another 
form of dark matter is assumed to exist, these models have 
been included, but their cross sections have been scaled by 
a factor fl^l ^CDm- In Fig. [8^, S^^m = 45 MeV has been 
chosen, and in[8j3 results for Z'ttAt = 64 MeV are shown for 
comparison. As one can see, the cross sections shift higher 
by a factor of a few for the larger value of St^n- 

Also shown on the plot are the current CDMS [56], 
XENONIO [57] and the newest CDMS [60] exclusion curves 
which place an upper limit on the scattering cross section. 
As one can see, the current limits have only just now be- 
gun to probe CMSSM models. 




10-" 

10-' 
10-10 

lo-'H 



10-1 



CMSSM, ^>0, 5:=45 MeV 



CDMS II 



CDMS II (2008) 








100 600 

(GeV) 



800 



10-= 



a. 



CMSSM, n>0, £=64 MeV 
CDMS II 




my, (GeV) 



Fig. 8. Scatter plots of the spin-independent elastic- scattering 
cross section predicted in the CMSSM for (a) X'vrjv = 45 MeV 
and (b) 64 MeV. 
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5 Indirect sensitivities 

Measurements of electroweak precision observables (EWPO) 
as well as B-physics observables (BPO) can provide in- 
teresting indirect information about the supersymmetric 
parameter space. We have already seen the impact of mea- 
surements of the anomolous magnetic moment of the muon, 
the branching ratio oi b sj in addition to the non- 
discovery of charginos and the Higgs boson at LEP which 
impose significant lower bounds on mi/2- 

The following EWPO were considered in [6T 1 l62 p 63j: 
the W boson mass, Mw, the effective weak mixing angle 
at the Z boson resonance, sin^ ^^ofr, the width of the Z, Fz, 
the Higgs mass, m/i, and (g^ — 2) in addtion to the BPO: 
h decays h ^ s^, ^ IJ'^ i Bu — > tVt and the mass 
mixing parameter AMb^- An analysis was performed of 
the sensitivity to mi/2 moving along the WMAP strips 
with fixed values of Aq and tan/3. The experimental cen- 
tral values, the present experimental errors and theoretical 
uncertainties are as described in [63]. Assuming that the 
nine observables listed above are uncorrelated, a fit has 
been performed with 



E 



R 



>thco 



2 log 



(32) 

Here R'!^^^ denotes the experimental central value of the 
nth observable {Mw^ sin^6'cff, -Tz, (ff — '2-)^ and BR(& 
57), BR(B„ TVr), AMb,), is the corresponding 

MSSM prediction and cr„ denotes the combined error. Ad- 
ditionally, cr™'" is the minimum combined error over the 
parameter space of each data set, and Xm^ ^^^^ Xb 
note the contribution coming from the experimental 
limits on the lightest MSSM Higgs boson mass and on 
BR(i?s — > respectively (63]. 

As examples of contributions to x^, predictions of some 
of the electroweak and B-physics observables are shown 
in Figures [H - [Til taken from [63J . The current status of 
the MSSM prediction and the experimental resolution is 
shown in each figure. For each value of j4o/mi/2, we plot 
two sets of points. For each mi/2, is chosen to lie in the 
co-anniliation strip (low values of mo) or the focus point 
region (high values of mg). 

The CMSSM predictions for My/ in the coannihilation 
and focus-point regions shown in Fig. [9] are quite similar, 
and depend little on Aq. We also see that small values 
of mi/2 are slightly preferred, refiecting the familiar fact 
that the experimental value of Mw is currently somewhat 
higher than the SM prediction. 

The corresponding results for sin^ 0off in the CMSSM 
are shown in Fig. [10] for tan/? = 10 as functions of mi/2- 
Best agreement occurs at large mi/2 values. However, tak- 
ing all uncertainties into account, the deviation for mi/2 
generally stays below the level of one sigma. We note that 
the predictions for sin^ O^a in the coannihilation and focus- 
point regions are somewhat different. 

The predictions for mn in the CMSSM for tan^ = 10 
are shown in Fig. [TT] The predicted values of m/i are sim- 
ilar in the coannihilation and focus-point regions. They 



CMSSM, n>0, 1 
tan(5 = 1 0, A(, = 

= 10,A„ = +m„2 
= 10,A„ = -m„^ 
! 10, A„ = +2 m„ 



tanp = 
tanp = 
tanp = 
tanp = 



10, A =-2 m„ 




m^,2 [GeV] 

Fig. 9. The CMSSM predictions for Mw are shown as func- 
tions of mi/2 along the WMAP strips for tan (3 = 10 for various 
A(i values. The centre (solid) line is the present central experi- 
mental value, and the (solid) outer lines show the current ±l-cr 
range. The dashed lines correspond to the full error including 
also parametric and intrinsic uncertainties. 



CMSSM, n > 0, rTlj = 1 71 .4 GeV 
tanp = 10, A;, = 
tanp = 10, Aj = +m,,2 
tanp = 10, Aj = -m,,2 
tanp = 10, Aj, = +2 m^^j, 
tanp = 10, Aj = -2 m,,2 



200 400 600 800 1000 

m„2 [GeV] 

Fig. 10. The CMSSM predictions for sin^ O^ff as functions of 
mi/2 along the WMAP strips for tan /? = 10 for various Ao val- 
ues. The centre (solid) line is the present central experimental 
value, and the (solid) outer lines show the current ±1-(t range. 
The dashed lines correspond to the full error including also 
parametric and intrinsic uncertainties. 



depend significantly on Aq, particularly in the coannihi- 
lation region, where negative values of Aq tend to predict 
very low values of that are disfavoured by the LEP 
direct search. Also shown in Fig. [Tl] is the present nomi- 
nal 95 % C.L. exclusion limit for a SM-like Higgs boson, 
namely 114.4 GeV [34J, and a hypothetical LHC measure- 
ment of ruh = 116.4 ± 0.2 GeV. 

The current status of the CMSSM prediction and the 
experimental resolution for the anomalous magnetic mo- 
ment of the muon, Aa^^ is shown in Fig. 1121 where the 1- 
and 2-cr bands are shown. We note that the coannihilation 
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Fig. 11. The CMSSM predictions for nih as functions of mi/2 
with tan/3 = 10 for various Aq. We also show the present 
95% C.L. exclusion limit of 114.4 GeV and a hypothetical LHC 
measurement of = 116.4 ± 0.2 GeV. 
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CMSSM,|j>0, mj = 171.4 GeV 
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tanp = 10, = +m^/2 
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tanp = 10, A„ = +2 

tanp = 10, A„ = -2 



Fig. 13. The CMSSM predictions for BR{b 37) as functions 
of mi/2 along the WMAP strips for tan/3 — 10 and various 
choices of Aq. The central (solid) line indicates the current 
experimental central value, and the other solid lines show the 
current ±l-cr experimental range. The dashed line is the il-a 
error including also the full intrinsic error. 



and focus-point region predictions for are quite differ- 
ent. For tan/3 = 10, the focus-point prediction agrees less 
well with the data. 
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Fig. 12. The CMSSM predictions for {g - 2)^, Aa^, as func- 
tions of mi/2 along the WMAP strips for tan /3 = 10 for various 
Ao values. In each panel, the centre (solid) line is the present 
central experimental value, and the solid (dotted) outer lines 
show the current ±l(2)-(j ranges. 



Next we turn to some B-physics observables. In Fig. [T3l 
we show the predictions in the CMSSM for BR{b — > 57) 
for tan/3 = 10 as functions of mi/2, compared with the 1- 
a experimental error (full line) and the full error (dashed 
line, but assuming a negligible parametric error). We see 
that positive values of Aq are disfavoured at small mi/2. 

Finally, in Fig.[Tlthe CMSSM predictions for BR{Bs 
fi'^fi~) for tan/3 = 10 as functions of toi/2 are compared 



with the present Tevatron limit. For tan /3 10, the CMSSM 
prediction is significantly below the present and future 
Tevatron sensitivity. However, already with the current 
sensitivity, the Tevatron starts to probe the CMSSM coan- 
nihilation region for tan (3 — 50 (not shown) and Aq > 0, 
whereas the CMSSM prediction in the focus-point region 
is significantly below the current sensitivity. 
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Fig. 14. The CMSSM predictions for BR{Bs as 
functions of mi/2 along the WMAP strips for tan/3 — 10 and 
various choices of Aq. The solid line indicates the current ex- 
perimental 95% C.L. exclusion bound. 



In Fig. [15] the net result of the analysis in the 
CMSSM is shown using the combined values for the 
EWPO and BPO, computed from Eq. ^ We see that 
the global minimum of ~ 4.5 for both values of tan /?. 
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This is quite a good fit for the number of experimental ob- 
servables being fitted. For both values of tan/3, the focus- 
point region is disfavoured by comparison with the coan- 
nihilation region, though this effect is less important for 
tan/3 = 50. For tan/3 = 10, mi/2 - 300 GeV and Aq > 
are preferred, whereas, for tan/3 = 50, mi/2 ~ 600 GeV 
and Aq < are preferred. This change-over is largely due 
to the impact of the LEP nih constraint for tan/3 = 10 
and the 6 — > 57 constraint for tan f3 ~ 50. 
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CMSSM, n > 0, rTij = 1 71 .4 GeV- 
tanp = 10, A„ = 
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CMSSM, n > 0, rrij 
tanp = 50, A„ = 
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tanp = 50, A„ = -2 m,. 



171.4 GeV- 



200 400 600 800 1000 1200 1400 1600 1800 2000 

m„2 [GeV] 

Fig. 15. The combined function for the electroweak ob- 
servables Mw, sin^^efi, ^z, {g — 2)fi, m/i, and the b physics 
observables BR(6 57), BR{Bs /i^M"), BR(B„ -* ru-r) 
and AMb,, evaluated in the CMSSM for tan/3 = 10 (a) 
and tan/3 = 50 (b) for various discrete values of Aq. We use 
mt = 171.4 ± 2.1 GeV and mb(mb) = 4.25 ± 0.11 GeV, and 
Too is chosen to yield the central value of the cold dark matter 
density indicated by WMAP and other observations for the 
central values of mt and mb{mb). 



6 GUT-less models 

The input scale at which universality is assumed in CMSSM 
models is usually taken to be the SUSY GUT scale, Mqut ^ 



2 X 10^^ GeV. However, it may be more appropriate in 
some models to assume the soft SUSY-breaking param- 
eters to be universal at some different input scale, Mi„. 
Specific scenarios in which the soft SUSY-breaking pa- 
rameters may be universal at a scale below Mqjjt occur 
in models with mixed modulus-anomaly mediated SUSY 
breaking, also called mirage-mediation 64], and models 
with warped extra dimensions [65^ . In the case of mirage- 
mediation, the universality scale is the mirage messenger 
scale, which is predicted to be Mi„ ~ 10^° — 10^^ GeV in 
the case of KKLT-style moduli stabilization . In other 
models, the universality scale may lie anywhere between 
1 TeV and Mpi. 

In the CMSSM with universality imposed at the GUT 
scale, the one-loop renormalizations of the gaugino masses 
Ma, where a = 1,2,3, are the same as those for the corre- 
sponding gauge coupHngs, Ua- Thus, at the one-loop level 
the gaugino masses at any scale Q < Mqjjt is given in 
Eq. [15] with t = Q. On the other hand, in a GUT-less 
CMSSM, where the gauge-coupling strengths run at all 
scales below the GUT scale but the soft SUSY-breaking 
parameters run only below the lower universality scale, 
Mi„, at which all the gaugino masses are assumed to be 
equal to mi/2 = Ma{Min), we have 



MaiQ) 



aa{Q) 



-mi/2 



(33) 



at the one-loop level. Since the runnings of the coupling 
strengths in GUT and GUT-less CMSSM scenarios are 
identical, the low-energy effective soft gaugino masses, 
Ma{Q), in GUT-less cases are less separated and closer 
to mi/2 than in the usual GUT CMSSM. This is seen in 
Fig. [16] where the low energy gaugino masses are shown 
as a function of the input unification scale. 

Similarly, the sfermion masses also run less and have 
masses closer to mo at the weak scale than their corre- 
sponding masses in the CMSSM. The soft SUSY-breaking 
scalar masses are renormalized by both gauge and (partic- 
ularly in the cases of third-generation sfermions) Yukawa 
interactions, so the running is somewhat more compli- 
cated. At the one-loop level one can summarize the effects 
of renormalizations at any Q < Min as 

ml^iQ) ^ ml{M,n) + C,(g, M,„)m2/2, (34) 

where we make the CMSSM assumption that the r7i§(Mi„) 
are universal at Mj„, and the Ci{Q, Min) are renormaliza- 
tion coefficients that vanish as Q ^ Min- The two-loop- 
renormalized soft SUSY-breaking masses of the the first- 
and second-generation left- and right-handed squarks, qL,R-i 
the stop mass eigenstates, ii,2, and the left- and right- 
handed sleptons, iL,R are shown in Fig.[T7] We see again 
that in GUT-less cases the soft SUSY-breaking scalar masses 
are less separated and closer to mo than in the usual GUT- 
scale CMSSM. 

One of the most dramatic changes when Mi„ is low- 
ered from the GUT scale is its effect on the footprint in 
the (m]^/2i"^o) plane of the constraint on the relic abun- 
dance of neutralinos inferred from WMAP. In the GUT- 
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Fig. 16. The dependences of observable gaugino mass param- 
eters on the input scale Mi„ at which they are assumed to be 
universal. The calculations are made for the representative case 
mi/2 = 700 GeV, mo = 1000 GeV, Ao = 0, tan/3 = 10 and 
fi>0. 
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Fig. 17. As in Fig. \16l the dependences of observable scalar 
mass parameters on the input scale Mi„ at which they are as- 
sumed to be universal. 



less CMSSM scenario [28], as the universality scale is low- 
ered to Alin ^ 10^^ GeV, the co-annihilation strip, the 
funnel and the focus point regions of the CMSSM ap- 
proach each other and merge, forming a small WMAP- 
preferred island in a sea of parameter space where the 
neutralino relic density is too small to provide all the cold 
dark matter wanted by WMAP. 

For comparison purposes, the CMSSM plane for tan /3 = 
10 is shown in Fig. [TH] There are already changes in the 
relic density as the universality scale is lowered to Mm = 



10^^ GeV seen in Fig. \T9\ The allowed focus-point region 
starts to separate from the LEP chargino bound, moving 
to larger mi/2- For Min — 10^^ GeV, the allowed focus- 
point region also dips further down, away from the elec- 
troweak vacuum condition boundary, while the coannihi- 
lation strip moves up and farther away from the region 
where the stau is the LSP. This is shown in Fig. [201 An- 
other remarkable feature at this value of Min is the ap- 
pearance of the rapid-annihilation funnel, familiar in the 
GUT-scale CMSSM at large tan /3, but an unfamiliar fea- 
ture for tan/3 = 10. 
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tan p = 10 , > 



> 
O 




1000 2000 

mi/2 (GeV) 

Fig. 18. The (mi/2,'Tio) plane with tan/3 = 10 and Ao = 
for Alin ~ Mgut ■ We show contours representing the LEP 
lower limits on the chargino mass (black dashed line), and the 
Higgs bound (red dot-dashed). We also show the region ruled 
out because the LSP would be charged (dark red shading), and 
that excluded by the electroweak vacuum condition (dark pink 
shading). The region favoured by WMAP has light turquoise 
shading, and the region suggested 67/(7^ — 2 at 2-a has medium 
(pink) shading, with the 1-a contours shown as black dashed 
lines. The solid (dashed) dark blue contours correspond to the 
approximate sparticle reach with 10 (TO) fb~^ of integrated 
LHC luminosity, as discussed in the text. 



As the universality scale is further decreased to Min = 
-|^q12.5 QgY^ a,s shown in panel Fig. [211 the atoll formed 
by the conjunction of what had been the focus-point and 
coannihilation strips has shrunk, so that it lies entirely 
within the range of (?tii/2, tiq) shown in the figure. We 
now see clearly two distinct regions of the plane excluded 
due to an excess relic density of neutralinos; the area en- 
closed by the atoll and the slice between the lower funnel 
wall and the boundary of the already-excluded r-LSP re- 
gion. 

In Fig. [22] for M„ = lO^^ GeV, the focus-point and 
coannihilation regions are fully combined and the atoll has 
mostly filled in to become a small island of acceptable relic 
density. To the right of this island is a strip that is pro- 
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Fig. 19. As m Fig.\Mfor Mi-a = lO" GeV. 
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Fig. 21. As in Fig.\Mfor Mi-a = 10^^'^ GeV. 
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Fig. 20. As m Fig.[Mfor Mi-a = lO" GeV. 
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Fig. 22. As m Fig. {MM Min = 10^^ GeV. 



vided by the lower funnel wall. The strip curves slightly 
as mi/2 increases then takes a sharp plunge back down 
towards the boundary of the region where the stau is the 
LSP, a feature associated with the XX h-\- A threshold. 
Reduction in the universality scale from this point results 
in the lower funnel wall being pushed down into the ex- 
cluded T LSP region and total evaporation of the island as 
seen in Figs . and [Ml and only a small residual turquoise 
region at large mi/2 where the relic density is within the 
WMAP limits. At all other points in the visible part of 
the (mi/2j Too) plane the relic density of neutralinos is too 
low to provide fully the cold dark matter density preferred 
by WMAP et al. Of course, these SUSY models would not 
be excluded if there is another source of cold dark matter 
in the universe. 



The discovery potential of ATLAS was examined in [57] , 
and more recently a CMS analysis ^68] has provided reach 
contours in the (mo, 7711/2) plane within the CMSSM for 
tan/3 = 10. Both studies found that the greatest discov- 
ery potential is achieved by an inclusive analysis of the 
channel with missing transverse energy, E^^'"^ , and three 
or more jets. To a good approximation, the contours de- 
pend only on m~ and to~, although processes involving 
other gauginos and sleptons may become important near 
the focus-point and coannihilation strips |69j . The 5-cr in- 
clusive supersymmetry discovery contours in the CMSSM 
for 1.0 and 10 fb~^ of integrated LHC luminosity are 
shown for tan/3 = 10 in Figure 13.5 of Ref. [68]. Since 
the inclusive reach is expected to be fairly linear above 
Too — 1-5 TeV, these contours have been extended [70] 
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Fig. 23. As in Fig.\Mfor M^^ = lO"'^ GeV. 
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Fig. 24. As m Fig.UEfor Mi-a = 10^" GeV. 

linearly above mo = 1200 GeV, then the sensitivity was 
fit with a third-order polynomial in mo and mi/2 to ex- 
tend the approximate LHC supersymmetry reach out to 
mo = 2 TeV, as shown in Fig. [25l These fits are compared 
with the CMS reaches. The largest differences between the 
approximate reach contours and the contours shown in the 
CMS TDR [68] are ~ 25 GeV for the 10 fb'^ contour and 
^ 50 GeV for the 1 fb~^ contour. 

The next step is to change variables from (mi/2, mo) 
{m~,m~) using ((33)) and (p4)) . Starting from the contours 
specified in Fig. [21] as functions of the gluino and squark 
masses, for each value of Min, the discovery contours were 
translated [7^ back into the corresponding (mi/2,mo) 
plane. Clearly, the contours move as the universality scale 
Min is lowered and the gluino and squark masses change 




200 400 600 800 1000 



mi/2 (GeV) 

Fig. 25. Approximate LHG supersymmetry reach contours 
for integrated luminosities of 1 fb~^ and 10 fh~^ (smooth 
curves), compared with the expected CMS reach given in the 
CMS TDR 168] for tan /? = 10 fyO[ /. 



according to ([55)1 and ([M)l . The shifted contours are dis- 
played appropriately in Figs. [18]- [2ll showing the change 
in the reach across the mi/2, mo plane with respect to the 
regions preferred cosmologically. 



7 mSUGRA models 

As discussed earlier, mSUGRA models, in contrast to the 
CMSSM, have an additional boundary condition, namely 
the value of the bilinear B-terni is fixed at the GUT scale 
to i?o = ^0 ~ 'Tio. As a consequence, one is no longer 
free to choose a fixed value of tan (3 across the mi/2 — mo 
plane. Instead, tan/3 is derived for each value of mi/2, mo 
and Aq [23 . Phenomenologically distinct planes are rather 
determined by a choice for Ao/rriQ. In Fig. [261 two such 
planes are shown with (a) Aq / mo = 3 — \/3 as predicted 
in the simplest model of supersymmetry breaking I27j and 
with (b) Ao/mo = 2.0. Shown are the contours of tan/? 
(solid blue lines) in the (m]^/2i"^o) planes. Also shown 
are the contours where m^± > 104 GeV (near-vertical 
black dashed lines) and ruh > 114 GeV (diagonal red 
dash-dotted lines). The regions excluded by & ^ 57 have 
medium (green) shading, and those where the relic den- 
sity of neutralinos lies within the WMAP range have light 
(turquoise) shading. The region suggested by — 2 at 2-a 
has very light (pink) shading. As one can see, relatively 
low values of tan /3 are obtained across the planes. 

Another difference between the CMSSM and models 
based on mSUGRA concerns the mass of the gravitino. 
In the CMSSM, it is not specified and and can be taken 
suitably large so that the neutralino is the LSP (outside 
the stau LSP region). In mSUGRA, the scalar masses at 
the GUT scale, mo, are determined by (and equal to) the 
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Fig. 26. Examples of mSUGRA {rni/2,mo) planes with con- 
tours of tan/3 superposed, for /i > and (a) the simplest 
Polonyi model with Ao/mo = 3 — -\/3, and (h) Ao/nio = 2.0, 
all with Bo = Ao — mo- In each panel, we show the regions 
excluded by the LEP lower limits on MSSM particles and those 
ruled out by b —* s'y decay (medium green shading): the regions 
favoured by — 2 are very light (yellow) shaded, bordered by 
a thin (black) line. The dark (chocolate) solid lines separate 
the gravitino LSP regions. The regions favoured by WMAP in 
the neutralmo LSP case have light (turquoise) shading. The 
dashed (pink) line corresponds to the maximum relic density 
for the gravitino LSP, and regions allowed by BBN constraint 
neglecting the effects of bound states on NSP decay are light 
(yellow) shaded. 



gravitino mass. In Fig. 1261 the gravitino LSP and the neu- 
tralino LSP regions are separated by dark (chocolate) sohd 
lines. Above this line, the neutralino (or stau) is the LSP, 
whilst below it, the gravitino is the LSP [7T]. As one can 
see by comparing the two panels, the potential for neu- 
tralino dark matter in mSUGRA models is dependent on 
Ao/mo. In panel (a), the only areas where the neutralino 
density are not too large occur where the Higgs mass is far 



too small or, at higher mo the chargino mass is too small. 
At larger Ao/mo, the co-annihilation strip rises above the 
neutralino-gravitino delineation. In panel (b), we see the 
familar co-annihilation strip. It should be noted that the 
focus point region is not realized in mSUGRA models as 
the value of ^jl does not decrease with increasing mo when 
Ao/mo is fixed and i?o = — "lo. There are also no 
funnel regions, as tan/3 is never sufficiently high. 

In the gravitino LSP regions, the next-to-lightest su- 
persymmetric particle (NSP) may be either the neutralino 
or stau which is now unstable. The relic density of grav- 
itinos is acceptably low only below the dashed (pink) line. 
This excludes a supplementary domain of the (mx/2,mo) 
plane in panel (a) which has a neutralino NSP (the dot- 
ted (red) curve in panel (a) separates the neutralino and 
stau NSP regions). However, the strongest constraint is 
provided by the effect of neutralino or stau decays on big 
bang neucleosynthesis (BBN) [7? . Outside the light (yel- 
low) shaded region, the decays spoil the success of BBN. 
Note that in panel (b), there remains a region which is 
excluded because the stau is the LSP. 

Recently, new attention has been focussed on the re- 
gions with a stau NSP due to its ability to form bound 
states (primarily with ■*He). When such bound states oc- 
cur, they catalyze certain nuclear reactions such as '*IIe(D, 
7)^Li which is normally highly suppressed due to the pro- 
duction of a low energy 7 whereas the bound state reaction 
is not [73]. In Fig. [27b. the (mi/2,mo) plane is displayed 
showing explicit element abundance contours [7J when 
the gravitino mass is m^/2 = 0.2mo in the absence of stau 
bound state effects. To the left of the solid black line the 
gravitino is the not the LSP. The diagonal red dotted line 
corresponds to the boundary between a neutralino and 
stau NSP. Above the line, the neutralino is the NSP, and 
below it, the NSP is the stau. Very close to this boundary, 
there is a diagonal brown solid line. Above this line, the 
relic density of gravitinos from NSP decay is too high, i.e.. 



"^3/2 

rriNSP 



^Nsph^ > 0.12. 



(35) 



Thus we should restrict our attention to the area below 
this line. 

The very thick green line labelled ^Li = 4.3 corre- 
sponds to the contour where ''Li/H = 4.3 x lO^^*', a value 
very close to the standard BBN result for ''Li/H. It forms 
a 'V shape, whose right edge runs along the neutralino- 
stau NSP border. Below the V, the abundance of ''Li is 
smaller than the standard BBN result. However, for rel- 
atively small values of mx/2 5 the ^Li abundance does not 
differ very much from this standard BBN result: it is only 
when mi/2 ^ 3000 GeV that ^Li begins to drop signifi- 
cantly. The stau lifetime drops with increasing mi/2, and 
when T ~ 1000 s, at mi/2 ~ 4000 GeV, the "^Li abundance 
has been reduced to an observation-friendly value close to 
2 X 10~^° as claimed in [75| and shown by the (unlabeled) 
thin dashed (green) contours. 

The region where the ^Li/^Li ratio lies between 0.01 
and 0.15 forms a band which moves from lower left to 
upper right. As one can see in the orange shading, there is 
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Fig. 27. The (mi/2,mo) planes for Aq — 0, fi > and 

tan/3 = 10 with 7713/2 = 0.2mo without (a) and with (b) the 
effects of boundstates included. The regions to the left of the 
solid black lines are not considered, since there the gravitino is 
not the LSP. In the orange (light) shaded regions, the differ- 
ences between the calculated and observed light-element abun- 
dances are no greater than in standard BBN without late par- 
ticle decays. In the pink (dark) shaded region in panel b, the 
abundances lie within the ranges favoured by observation. The 
significances of the other lines and contours are explained in 
the text. 



a large region where the lithium isotopic ratio can be made 
acceptable. However, if we restrict to D/H < 4.0 x 10~^, 
we see that this ratio is interesting only when ^Li is at or 
slightly below the standard BBN result. 

Turning now to Fig. 127b . we show the analogous re- 
sults when the bound-state effects are included in the cal- 
culation. The abundance contours are identical to those in 
panel (a) above the diagonal dotted line, where the NSP is 
a neutralino and bound states do not form. We also note 



that the bound state effects on D and ^He are quite mini- 
mal, so that these element abundances are very similar to 
those in Fig. [27k. However, comparing panels (a) and (b), 
one sees dramatic bound-state effects on the lithium abun- 
dances. Everywhere to the left of the solid blue line labeled 
0.15 is excluded. In the stau NSP region, this means that 
Wi/2 ^ 1500 GeV. Moreover, in the stau region to the 
right of the ^Li/^Li = 0.15 contour, the ^Li abundance 
drops below 9 x 10^^^ (as shown by the thin green dot- 
ted curve). In this case, not only do the bound-state effects 
increase the ^Li abundance when is small (i.e., at rel- 
atively long stau lifetimes), but they also decrease the ^Li 
abundance when the lifetime of the stau is about 1500 s. 
Thus, at (mi/2,mo) ~ (3200,400), we find that ^Li/'^Li 
~ 0.04, ■^Li/H ~ 1.2 X 10-1°, and D/H ~ 3.8 x 10^. Indeed, 
when mi/2 is between 3000-4000 GeV, the bound state ef- 
fects cut the ^Li abundance roughly in half. In the darker 
(pink) region, the lithium abundances match the observa- 
tional plateau values, with the properties ^Li/^Li > 0.01 
and 0.9 X lO"!" < ^Li/H < 2.0 x 10-i°. 

Next we return to an example of an mSUGRA model 
based on the Polonyi model for which /rriQ = 3 — 
with the condition that 7713/2 = "niQ- In Fig. 128b ,. we show 
the mSUGRA model without the bound states. In the 
upper part of the plane, we do not have gravitino dark 
matter. We see that '^He/D [76] eliminates all but a tri- 
angular area which extends up to mo = 1000 GeV, when 
mi/2 = 5000 GeV. Below the ^He/D = 1 contour, D and 
^Li are close to their standard BBN values, and there is 
a substantial orange shaded region. We note that ^Li is 
interestingly high, between 0.01 and 0.15 in much of this 
region. 

As seen in Fig. [28b . when bound-state effects are in- 
cluded in this mSUGRA model, both lithium isotope abun- 
dances are too large except in the extreme lower right 
corner, where there is a small region shaded orange. How- 
ever, there is no region where the lithium abundances fall 
within the favoured plateau ranges. 

The prospects for direct detection in the mSUGRA 
models (when the neutralino is the LSP) is very similar to 
that of the CMSSM but shghtly more difficult. In a plot 
such as that shown in Fig. [H the mSUGRA predictions 
are down by a factor of about 2 or less. 



8 The NUHM 

In the NUHM, the Higgs soft masses are treated inde- 
pendently from mo. In effect, this allows one to choose 
^ and the Higgs pseudoscalar mass, niA freely (up to 
phenomenological constraints). Two examples of planes 
in the NUHM are shown in Figs. [M] and [301 In Fig. 
an mi/2 — mo plane is shown with /x = 700 GeV and 
ruA — 400 GeV fixed across the plane [30]. As usual, the 
light (turquoise) shaded area is the cosmologically pre- 
ferred region. There is a bulk region satisfying this prefer- 
ence at mi/2 - 150 GeV to 350 GeV and mo - 100 GeV. 
The dark (red) shaded regions are excluded because a 
charged sparticle is lighter than the neutralino. As in the 
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Fig. 28. As in Fig.\2^ for the mSUGRA case with Ao/mo 
3-^3. 



CMSSM, there are light (turquoise) shaded strips close to 
these forbidden regions where coannihilation suppresses 
the relic density sufficiently to be cosmologically accept- 
able. Further away from these regions, the relic density is 
generally too high. At small mi/2 and mo the left handed 
sleptons, and also the sneutrinos, become lighter than the 
neutralino. The darker (dark blue) shaded area is where 
a sneutrino is the LSP. 

The near-vertical dark (black) dashed and light (red) 
dot-dashed lines in Fig.[29]are the LEP exclusion contours 
m^± > 104 GeV and rrih > 114 GeV respectively. As in 
the CMSSM case, they exclude low values of and 
hence rule out rapid relic annihilation via direct-channel 
h and poles. The solid lines curved around small val- 
ues of mi/2 and mp bound the light (pink) shaded region 
favoured by and recent analyses of the e+e^ data. 



strip with low fi^. 
density on either side. The low-f?-^ 



A striking feature in Fig. [29] when mi/2 ^ 450 GeV is a 
h"^ , which has bands with acceptable relic 
strip is due to rapid 
annihilation via the direct-channel A, H poles which occur 
when — mA/'2. = 200 GeV, indicated by the near- 
vertical solid (blue) line. These correspond to the funnel 
regions found in the CMSSM [l9pT]. but at larger tan/3. 
At higher mi/2 ~ 1300— 1400 GeV, there is another region 
of acceptable relic density. This band, the transition band, 
is broadened because the neutralino acquires significant 
Higgsino content as r7T,i/2 becomes greater than /x, and 
the relic density is suppressed by the increased W^W~ 
production. To the right of this band, the relic density 
falls below the WMAP value. 



tan P = 10 , n = 700 GeV, = 400 GeV 



> 




mi/2 (GeV) 

Fig. 29. Projections of the NUHM model on the (mi/2,'Tio) 
planes for tan/? = 10 and n — 700 GeV and rriA = 400 GeV. 
The (red) dot-dashed lines are the contours mn = 114 GeV as 
calculated using FeynHiggs \35^ . and the near-vertical (black) 
dashed lines are the contours m^± = 104 GeV. The light 
(turquoise) shaded areas are the cosmologically preferred re- 
gions. The dark (brick red) shaded regions is excluded because 
a charged particle is lighter than the neutralino, and the darker 
( dark blue ) shaded regions is excluded because the LSP is a 
sneutrino. There is a very small medium (green) shaded region 
excluded by b —> 37, at small mi/2. The regions allowed by the 
E821 measurement of at the 2-a level, are shaded (pink) 
and bounded by solid black lines. 



Another example of an NUHM plane is shown in Fig. 
1301 representing a /i — ra^ plane with a fixed choice of 
values of tan/3 = 35, mi/2 = 500 GeV, toq = 1000 GeV 
and = [58]. Line types and shading are as in Fig. 
[29l The region with acceptable relic density now takes the 
form of a narrow strip of values of /i ^ ±300 350 GeV 
where the relic density lies within the WMAP range for 
almost all values of m^. The exception is a narrow strip 
centred on tua ~ 430 GeV, namely the rapid-annihilation 
funnel where ~ mA/'i-, which would be acceptable if 
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there is some other source of cold dark matter. Note the 
increased importance of the — > 57 constraint which now 
excludes almost all of < 0. 

tan p = 35 , my2 = 500 , mo = 1000 , Aq = 




2000 



ii (GeV) 

Fig. 30. As in \2!A for the fi — tua plane with mi/2 = 500 
GeV, mo — 1000 GeV, tan/3 — 35 and Ao = 0. Here note that 
the constraint from b ^ sj excludes most of ^ < 0. 



An interesting representation of the NUHM parameter 
space is found in the construction of a (M^,tan/3) plane. 
A generic {Ma, tan (3) plane for fixed = 600 GeV, 

mo = 800 GeV, = 1000 GeV and Aq = [771i78j is shown 
in Fig. Ell The relic LSP density satisfies the WMAP con- 
straint only in narrow, near-vertical (pale blue) shaded 
strips crossing the plane. These lie to either side of the 
vertical (purple) line where = Myi/2. Within the nar- 
row unshaded strip straddling this line, the relic density 
is suppressed by rapid direct-channel annihilations to a 
value below the lower limit of the range for the cold dark 
matter density indicated by WMAP et al. This entire strip 
would be acceptable for cosmology if there were some addi- 
tional component of cold dark matter. Outside the shaded 
WMAP-compatible strips, at both larger and smaller val- 
ues of Ma, the relic LSP density is too high, and these 
regions are unacceptable. The LEP lower limit on ruh ex- 
cludes a strip of this plane at low Ma and/or tan/? indi- 
cated by the dash-dotted (red) line. In addition, (pink 
shading) prefers relatively large tan/3 > 36, while 6 — > 57 
excludes a (green shaded) region at low Ma and tan /?. 

It is clear from this example and the properties found 
in the /i — niA plane discussed above, that one can con- 
struct a TTiA — tan/? plane such that (nearly) the entire 
plane respects the cosmological relic density constraint 
[631178]. For example, from Fig. [3T1 we see that if one 
adjusts TOi/2 continuously as a function of Ma so as to 
remain within one of the narrow WMAP strips as Ma 
increases we can cover the entire rriA — tan /3 plane. Al- 
ternatively, we can use Fig. [301 as a guide to construct 
tua — tan (3 planes as well. Fixing /i in the range of 250 - 
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Fig. 31. A (MAytwP) plane with = 600GeV, mo — 

800 GeV, = 1000 GeV and Ao = 0. The range of cold dark 
matter density preferred by WMAP and other observations is 
attained in two narrow (pale blue) strips, one on either side 
of the vertical solid (blue) line where — Ma/2- The dark 
(green) shaded region at low Ma and low tan 13 is excluded by 
b — > S7, and the medium (pink) shaded region at tan/3 > 36 is 
favoured by a^. The region below the (red) dot-dashed line is 
excluded by the LEP bounds on nih- 



400 GeV will allow us to obtain a good relic density over 
the entire m^i — tan /? plane expect for the strip centered on 

= Ma/2. These planes can be used as particularly use- 
ful benchmarks for existing and future searchs for Higgs 
bosons [78]. We will consider two such planes which will 
be referred to as planes PI and P3. 

The plane PI is defined by fixed /i = 1000 GeV and 
mo = 800 GeV, while P3 is defined by fixed mi/2 = 500 
GeV and mo = 1000 GeV. As one can see, the light shaded 
(turquoise) region now fills up each plane except for a 
vertical swath in P3 where s-channel annihilation through 
rriA drives the relic density to small values. A likelihood 
analysis of these NUHM benchmark surfaces, including 
the EWPO Mw, sin^6'off, Fz, {g - 2)^ and m^ and the 
BPO BR(6 S7), BR(Bs ^ M'^M"), BR(B„ rvr) and 
AMb^ was performed recently in |63J . The lowest value 
in plane PI (P3) is x^j^ = 7.1(7.4) and the corresponding 
best-fit point is shown by a (red) cross in Fig.[321a'S are the 
= 2.30 and 4.61 contours around the best-fit point. 
The (black shaded) region in the lower left is excluded at 
the 95 % C.L. by the LEP Higgs searches in the channel 
e+e- ^ Z* ^ Zh,H 

The planes shown in Fig. [32] also show the 5-cr discov- 
ery contours for bb H/A — > r+r" at the LHC, where 
the r's decay to jets and electrons or muons. The analysis 
is based on 60 fb~^ for the final state r+r" —f jets [79] 
and on 30 fb^^ for the r+r" ^ e-l- jet [501 and r+r" 
/J, -I- jet [HI] channels, collected with the CMS detector. As 
shown in '821, the impact of the supersymmetric param- 
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Fig. 32. The (MA,tan;9) planes for the NUHM benchmark 
surfaces PI and PS, displaying the contours of Ay^ found in 
a recent global fit to EWPO and BPO f63f . All surfaces have 
Ao = 0. The dark shaded (black) region corresponds to the 
parameter region that is excluded by the LEP Higgs searches. 
Also shown are the 5- a discovery contours for H/A — > r^r^ 
at the LHC with 60 or 30 fb~^ (depending on the r decay chan- 
nels) and for t^u detection in the CMS detector when 



eters other than Ma and tan (3 on the discovery contours 
is relatively small in this channel. We see that the whole 
Zix^ < 2.30 regions of the surface PI would be covered by 
the LHC H/A — > r+r^ searches, and most of the corre- 
sponding regions of the surface P3. 

We also show in Fig. [321 the S-cr contours for discovery 
of the via its t^v decay mode at the LHC, in the 
case Mh± > mt. We see that the coverage is limited in 
each of the scenarios PI and P3 to Ma < 300 GeV and 
tan/3 > 30, reaching a small part of the Ax^ < 2.30 region 



of surface P3, and only a small part of the Ax^ < 4.61 
region of surface PI. One may also search for ^ t^u 
for lighter Mh± < rrit, but in the cases of surfaces PI this 
would be useful only in the regions already excluded by 
LEP, and the accessible regions in surfaces P3 would also 
be quite limited. 

In Fig. [33] the current reach of the XENONIO experi- 
ment fSTj is shown on planes PI and P3. The solid line in 
the figure correspond to the XENONIO bound obtained 
assuming St^n = 45 MeV. The corresponding reach ob- 
tained if S^N = 64 MeV is shown by the dotted curve. In 
plane PI, the current reach is minimal and does not enter 
the Ax^ < 4.61 region, whereas for plane P3 the best fit 
point would already have been probed by the XENONIO 
experiment if Z'^at — 64 MeV. If the strangeness content 
of the proton were zero, it would correspond to S-^n = 36 
MeV, and the XENONIO reach in this case is displayed 
by the dashed curve. Finally, we note that a future sensi- 
tivity to a cross section of 10~® pb would cover the entire 
surface in P3 while in PI would only extend slightly past 
the Z'^AT = 64 MeV curve. 

As can be inferred from the benchmark surface shown 
above, in the NUHM, current constraints already exclude 
many interesting models |58j . These models generically 
have large scattering cross sections. Scatter plots of show- 
ing the reach of CDMS and XENONIO for both E^n = 45 
MeV and 64 MeV are shown in Fig. [34] In comparison with 
Fig. [51 one sees that many more viable models are already 
currently being probed. 



9 A Hint of Higgs 

The CDF and DO Collaborations have already established 
important limits on the (heavier) MSSM Higgs bosons, 
particularly at large tan/3 [55 1 [51 1 iF5 1 l5 ^ [57] . Recently the 
CDF Collaboration, investigating the channel 

pp^0^r+T-, {<P^h,H,A), (36) 

has been unable to improve these limits to the extent of 
the sensitivity expected with the analyzed integrated lu- 
minosity of 1 fb~^ [86], whereas there is no indication 
of any similar effect in DO data [H]. Within the MSSM the 
channel p6p is enhanced as compared to the correspond- 
ing SM process by roughly a factor of tan^ /3 /((I + Ap)'^ + 
9) [88], where Af, includes loop corrections to the (pbb ver- 
tex (see [55] for details) and is subdominant for the r+r" 
final state. Correspondingly, the unexpected weakness of 
the CDF exclusion might be explicable within the MSSM 
if ruA w 160 GeV and tan /3 > 45. 

In the CMSSM, there are no parameter sets consistent 
with all the experimental and theoretical constraints. In 
particular, low values of toi/2 are required to obtain low 
values of m^. For tua — 160 GeV and the other parame- 
ters tuned to fulfill the B physics constraints, there are no 
CMSSM solutions for which the lightest Higgs mass is not 
in direct contradiction with experimental limits [34] . How- 
ever, we do find that all the constraints may be satisfied 
in the NUHM [89l[90]. 
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Fig. 33. The same {Ma, tan/3) planes for the NUHM bench- 
mark surfaces (a) PI and (h) P3 as in Fig. \3Sl displaying the 
expected sensitivities of present direct searches for the scatter- 
ing of dark matter particles. 



We consider first the {mi/2, ^o) plane shown in panel 
(a) of Fig. Eg which has tan/3 = 45, = 370 GeV and 
^0 = -1800 GeV, as well as Ma = 160 GeV. As be- 
fore, the dark (brown) shaded region at low toq is for- 
bidden, because there the LSP would be the lighter stau. 
The WMAP cold dark matter constraint is satisfied only 
within the lighter (turquoise) shaded region. To the left 
of this region, the relic density is too small, due to s- 
channel annihilation through the Higgs pseudoscalar A. 
As mi/2 increases away from the pole, the relic density 
increases toward the WMAP range. However, as rni/2 is 
increased, the neutralino acquires a larger Higgsino com- 
ponent and annihilations to pairs of W and Z bosons 
become enhanced. To the right of this transition region, 
the relic density again lies below the WMAP preferred 
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Fig. 34. Scatter plots of the spin-independent elastic- 
scattering cross section predicted in the NUHM (a) for St^n = 
45 MeV (b) and 64 MeV. 



value. The shaded region here is therefore an overlap of 
the funnel and transition regions discussed in |30j . The 
BK{Bs fi^) constraint [91] is satisfied between the 

outer black dash-dotted lines, labelled 10"'', representing 
the current limit on that branching ratio. Also shown are 
the contours where the branching ratio is 2 x 10~®, close 
to the sensitivity likely to be attainable soon by CDF and 
DO. Between these two contours, there is a strong can- 
cellation between the flavor-violating contributions aris- 
ing from the Higgs and chargino couplings at the one-loop 
level and the Wilson coefficient counterterms contributing 
to BR{Bs [77]. 

The dash-dotted (red) line shows the contour corre- 
sponding to nifi — 114 GeV, and only the region to the 
right of this line is compatible with the constraint im- 
posed by ruh (it should be kept in mind that there is 
still a ~ 3 GeV uncertainty in the prediction of [35]). 
Also shown in pink shading is the region favoured by 



(5 



at the two-(T level. The one- and two-cr con- 
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tours for {g — 2)^ are shown as elliptical dashed and solid 
black contours, respectively. The region which is com- 
patible with the WMAP relic density and rrih, and is 
also within the two-cr {g — 2)^ experimental bound, has 
BR(i?s fj,'^fj,~) > 2 X 10"^. The measured value of 
BR(6 — + 57) is in agreement with the theory prediction 
only to the left of the solid (green) region. 



lO-'i jilO-" /2.10-* 
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o 1000 
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300 1000 1500 
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Fig. 35. The NUHM parameter space as a function of mi/2 
and mo for fi = 370(400) GeV and Ao = -1800(-2000) GeV 
in panels a (b). We fix Ma = 160 GeV, tan/? = 45, mt = 
171.4 GeV and nib = 4.25 GeV. For the description of the var- 
ious lines and shaded areas, see the text. 



We see that there is a narrow wedge of allowed param- 
eter space in Fig. [35Fa). which has mi/2 ^ 600 GeV and 
mo 700 to 1100 GeV. The BR(6 sj) constraint is 
satisfied easily throughout this region, and {g — 2)^ cuts 
off the top of the wedge, which would otherwise have ex- 
tended to mo ^ 1500 GeV. Within the allowed wedge, 
mil is very close to the LEP lower limit, and BR(i?s — + 
) > 2 X 10"*. If Ma were much smaller (< 130 GeV), 



there would be no wedge consistent simultaneously with 
the ^^cDM, iTT-h and BR(i3s ^i^ ^T) constraints. 

Given the matrix element uncertainties for direct de- 
tection summarized above, we show in Fig. [35k ) the 1-cr 
lower limit on the calculated value of the elastic cross sec- 
tion as compared to the CDMS upper limit In the 
portion of the plane to the left of the (orange) dotted line, 
the lower limit on the calculated spin-independent elas- 
tic cross section is smaller than the CDMS upper bound, 
assuming the canonical local density. Whilst we have as- 
sumed St^n = 45 MeV, the calculated lower limit effec- 
tively assumes zero strangeness contribution to the proton 
mass, i.e., y = 0. In the region of interest, the lower limit 
on the calculated cross section is about 80% of the CDMS 
upper bound, whereas with a strangeness contribution of 
?; = 0.2, the cross section would exceed the CDMS bound 
by a factor of ~ 3. Thus, if Nature has picked this corner 
of the NUHM parameter space, we expect direct detec- 
tion of dark matter to be imminent. Consistency with the 
XENONIO limit would further require a reduction in the 
local dark matter density (to its lower limit). Intriguingly, 
the XENONIO experiment has seen some potential signal 
events that are, however, interpreted as background. 

Panel (a) of Fig. [551 explores the (/i, Ao) plane for the 
choice (mi/2,mo) = (600, 800) GeV, values close to the 
lower tip of the allowed wedge in Fig. [35ja). In this case, 
the region allowed by the BR(i?s /i+/j,~) constraint 
is helow the upper dash-dotted black line, and the LEP 
m/i constraint is satisfied only above the dash-dotted red 
line. We see that only a restricted range 360 GeV < fj, < 
390 GeV is compatible with the dark matter constraint. 
This corresponds to the transition strip where the neu- 
tralino is the appropriate bino/Higgsino combination. To 
the left of this strip, the relic density is too small and to 
the right, it is too large. Only a very restricted range of 
Aq ^ —1600 GeV is compatible simultaneously with the 
Tiih and BR(_Bs /i+/i^) constraints. Very large negative 
values of Ao are excluded as the LSP is the lighter stau. 
On the (/i, Aq) plane, the elastic scattering cross section is 
a rapidly decreasing function of /i and is almost indepen- 
dent of Aq. Indeed, NUHM points excluded by CDMS (or 
XENONIO) generally have low values of and Ma p5] . 
Values of /i > 355 GeV are compatible with CDMS if the 
strangeness contribution to the proton mass is negligible. 
For this choice of parameters, the entire displayed plane 
is compatible with BR(6 57) and [g — 2)^. 

Panel (b) of Fig. 1361 shows what happens if mx/2 is in- 
creased to 700 GeV, keeping toq and the other inputs the 
same. Compared to Fig. |36ja), we see that the WMAP 
strip becomes narrower and shifts to larger fji ~ 400 GeV, 
and that BR(5 sj) starts to exclude a region visible at 
smaller fi. If mi/2 were to be increased much further, the 
dark matter constraint and BR(6 sj) would no longer 
be compatible for this value of mo. We also see that, by 
comparison with [SSl^a). the BR{Bs IJ-'^IJ'~) constraint 
has moved to lower Aq, but the m/j constraint has dropped 
even further, and m/j > 114 GeV over the entire visible 
plane. The net result is a region compatible with all the 
constraints that extends from Aq ~ —1850 GeV down to 
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(to the left of the (orange) dotted hne) provided the strangeness 
contribution to the proton mass is smaU. 

Clearly our anticipation for the discovery of supersym- 
metry or new physics is great. If the CDF hint of the Higgs 
boson is realized, we should expect to see numerous de- 
partures from the standard model including the decay of 
Bg to n'^fJ.^ , BR(6 — > 57) should show deviations from the 
standard model value, and dark matter should be discov- 
ered by CDMS or XENONIO in the near future. Similarly 
if an observation of new physics in any one of these areas 
occurs, we should be ready for flood of new data showing 
departures from the standard model. 
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Fig. 36. The NUHM parameter space as a function of fj, and 
Aq for mi/2 = 600(700) GeV and mo = 800 GeV in panels a 
(b). We fix Ma = 160 GeV, tan/3 = 45, mt = 171.4 GeV and 
mb = 4.25 GeV. For the description of the various lines and 
shaded areas, see the text. 



Aq ^ —2150 GeV for fi ^ 400 GeV, with a coannihilation 
filament extending to larger /i when Aq ^ —2200 GeV. 
Once again all of the WMAP strip in this panel is com- 
patible with CDMS. 

The larger allowed area of parameter space is reflected 
in panel (b) of Fig. [35l which has = 400 GeV and = 
-2000 GeV, as well as tan/3 = 45 and Ma = 160 GeV as 
before. In this case, we see that a substantial region of the 
WMAP strip with mi/2 ^ 700 GcV and a width 5mi/2 ~ 
100 GeV, extending from mo ~ 750 GeV to higher mo is 
allowed by all the other constraints. The BR(Bs /i+/z~) 
and vrih constraints have now moved to relatively low val- 
ues of mi/2, but we stiU flnd BR(Bs — > > 2 x 10~* 



and rrih close to the LEP lower limit. The {g — 2)^ con- 
straint truncates the allowed region at toq ^ 1050 GeV. 
Once again, the allowed region is compatible with CDMS 
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